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Abstract 

Using high temperature expansions for the equal time correlator S{q) and 
static susceptibility x{q) for the t-J model, we present evidence for quantum 
critical (QC), z = l, behavior at intermediate temperatures in a broad range 
of t/J ratio, doping, and temperatures. We find that the dynamical sus- 
ceptibility is very close to the universal scaling function computable for the 
asymptotic QC regime, and that the dominant energy scale is temperature. 
Our results are in excellent agreement with measurements of the spin-echo de- 
cay rate, I/T2G, in La2Cu04, and provide qualitative understanding of both 
1/Ti and I/T2G nuclear relaxation rates in doped cuprates. 
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Recent interest in the doped antiferromagnets is related to the problem of high temper- 
ature superconductivity. While the magnetic behavior of the parent insulating compound, 
La2Cu04, can be described in terms of the S=l/2 Heisenberg model with the dominant 
interaction being the in-plane nearest neighbor exchange coupling J ~ 1500K, current un- 
derstanding of the doped materials is far from complete. The consensus on the details of 
the microscopic model which would quantitatively describe the magnetic properties over the 
entire doping range from the insulator to the fully doped compounds has yet to be reached. 

Recently, significant progress has been made in understanding the low energy spin dy- 
namics of these systems from a scaling and renormalization group point of view. As shown by 
Chakravarty, Halperin, and Nelson [|l|, the spin dynamics of an insulator, such as La2Cu04, 
is well described by the quantum nonlinear sigma (QNLcr) model. In case when the average 
sublattice magnetization is present at T = 0, the low temperature renormalized classical 
(RC) phase is characterized by an exponentially increasing antiferromagnetic correlation 
length, ^/a ~ exp(27rps/7'), where a is the lattice constant and ps spin stiffness (below we 
assume the units where h = ks = a = 1). Beyond the critical point, i.e., when the ground 
state does not possess Neel order, the quantum disordered (QD) phase has a finite ^ at T=0. 

Another aspect of the phase diagram of Ref. [|I|, which did not attract much attention 
until the recent work of Sachdev and Ye P], and Chubukov and Sachdev [§, is the quantum 
critical (QC) region, where in the leading order ,^~1/T Only at the critical point, i.e. 
at the boundary between the zero temperature Neel and disordered phases, QC behavior 
holds down to T = 0. In this case, ps vanishes, so that the model does not possess any energy 
scale, which thus is set by the temperature. For doped antiferromagnets, the possibihty of 
z=l criticality has been first pointed out in Ref. 0. 

Although for small temperatures a fine-tuning of the model parameters to the critical 
point is necessary to have the QC phase with C, ~ c/T, the range of parameters where it 
exists rapidly widens as the temperature increases (c is the spin wave velocity). While 
in the continuum QNLo" model the quantum critical behavior persists for arbitrarily high 
temperature, on the lattice the range of its applicability is necessarily limited form above, 
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i.e. quantum critical regime is intermediate between low and high temperatures. It has been 
argued in Ref. that the region of its applicability in the 2D S=l/2 Heisenberg model exists 
around T ~ 0.5J and that small doping should lead to a decrease in p^, thus extending the 
temperature range of QC behavior. 

Recently, it has been shown by Pines and one of us (A.S.) [^] using purely scaling con- 
siderations, that the experimental data of on the nuclear magnetic relaxation rates 
1/Ti and I/T2G in the superconducting cuprates imply a QC behavior at high tempera- 
tures over an unexpectedly broad doping range. This indicates that the high energy spin 
waves may survive even substantial doping, which would have important implications for 
superconductivity. 

In view of the above discussion, we present a study of the quantum-critical behavior 
in the doped 2D t-J model using the high temperature series expansion approach [^]. The 
nearest-neighbor version of the t-J model is described by the Hamiltonian: 

H = -tY. P{clc,^ + h.c.)P + J ^ S,S,-. (1) 

(P is the projection operator prohibiting double occupancy). The 10-term series in (3 = 
1/T for the equal time correlation function, ^(q), and static susceptibility, x(ct), has been 
generated for arbitrary t/ J, number of electrons per unit cell p = 1—x, as well as for arbitrary 
q. The substitution of the expansion variable w = tanh(/?//5o), which eliminates the influence 
of any singularities outside the stripe |Im/5| < 7r/5o/2, has been applied in order to improve 
the convergence; several different values of jSo have been used and the consistency of the 
approximations have been checked. Although our qualitative results hold for a broad range 
of t/J ratio, we chose to present our data for one particular value of t/J=l because the 
accuracy improves for smaller t/J. 

Our first important result is that in a broad range of the doping x, temperature, and t/J 
ratio, S{q) and x(q) near Q = (tt, tt) are well described by the following scaling expressions: 

= Sq ■ Sm{q^m), X(q) = Xq ■ Xm(gO> (2) 
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where q = (g^ universal scaling functions, while Sq, xq^ ^m, i'm 

depend on T, x, and t/J. We define Sm{i^) and Xm(/t) so that 5*^(0) = Xm(0) = 1 and 
dPSml dK^\f^=Q = (Px^/dn'^li^^o = —2. In this case, and correspond to the "second mo- 



ment" definition of the correlation length for 5'(q) and x(q)) respectively. Gradual character 
of deviation from the universal behavior described by Eq.(|D does not allow us to define an 
unambiguous boundary at which Eq.(0) fails. For t/J=l, the approximate range where the 
scaling behavior holds well is T=0.6J-J, x = — 15%; more detailed discussion will be given 
in a subsequent publication [|T^ . 

The scaling functions Sm and Xm are computable in the QNLo" model. We observe that 
the 1/N calculation of Ref. indicates negligible 1/N corrections to Xm and Sm in the QC 
regime (note that 1/N corrections to the prefactors and ^ are not negligible). Therefore, 
we can safely use scaling functions calculated for = oo, which we plot in Fig.|l] along 
with our data. We find that calculated 5*^ and Xm are very close to those obtained from 
the collapsed data for the t-J model (Fig.|lD. Further, in the asymptotic QC regime, 1/N 
calculation predicts that x{K)m is nearly Lorentzian, while S{K)m larger than Lorentzian for 
any given k. On the other hand, in the asymptotic RC regime, S{K,)m and x(^)m coincide 
and both are smaller than Lorentzian. We therefore conclude that the scaling functions of 
Eq.(^ indicate quantum critical behavior in the t-J model at moderate temperatures. 

We now turn to the dynamical properties. Since 



where g{C) = |^/ tanh (|C)) the ratio S'(q)/Tx(q) refiects the frequency distribution of the 




spectral weight at q. We find that for zero doping this ratio at q= (tt, vr) varies less than 1% 
in the range T=0.6J-J and is in excellent agreement with the 0(N) sigma model calculation 




(3) 




of Ref. 0: 
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in the asymptotic RC regime, this ratio would be equal to unity. In the doped case, the 
ratio somewhat increases, up to 1.16 for 15% doping at T = J. This increase may be caused 
by a broad electron-hole continuum, which does not significantly modify x"(q, cj) for small 
frequencies, but can yield significant contribution to the ratio because g{Q is large for ^ 1. 

Another universal quantity which may be temperature independent only when a) ~ T, 
i.e. in case of QC behavior, is 1 — Cm/ Cm- -^^^ insulator, this quantity indeed has less 
than 7% variation in the range T=0.6J-J. The comparison to 0{N) calculations of Ref. |^ 
for the asymptotic QC regime of the QNLcr model yields: 

I 0.043 ± 0.003 (t-J), 
^™ I 0.043 (N=3) and 0.035 (N=oo). 

In the doped case, Cm/ Cm becomes 10% smaller at T=0.6J, a decrease that may also be 
caused by the electron- hole continuum; it remains nearly doping independent at T=J. 

On the basis of the above arguments, we conclude that the dynamical susceptibility of 
the t-J model can be written as: 

X{^i^) =Xq-x(qC^ (6) 

where to the leading order uj T . In addition, in the asymptotic QC regime of the QNLa 
model, not only Eq.(|^) holds, but it is also expected that xq ~ C'^~^i where the critical 
exponent rj is nearly zero and can be neglected, and 1/^~T. We first calculate Xq/C"^ by 
generating 9-term series directly for this quantity and find that for T=0.6J-J it varies not 
more than 16%, compared to far greater (nearly by the factor of five in the undoped case) 
change in xq and separately. 

We now turn to the temperature dependence of C and evaluate it by generating 9-term 
series for ^T^/^ and 10-term series for Sq and XqT- C~^iT,x) for t/J =1 is plotted as a 
function of temperature in Fig.|^. As one can see, in a broad range of doping and temperature 
is nearly linear in T with doping independent slope. As shown by Chakravarty [|TT|, one 
would expect the temperature dependence to be of the form: 
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^ ~ ^ - C(x) T^-i/^ (7) 

with the critical exponent v ~ 0.7. The second term varies slowly as a function of tempera- 
ture over the interval of comparisons, and provides a doping dependent intercept for nearly 
linear ^^^(T). We note, however, that the slope of ^~^(T) determined from the numerical 
data differs nearly by the factor of two from the value calculated using Eq.(^ with T=0 
value of c and h ~ 0.8 — 1 determined in Refs. [p],p|,p!2[ . We conjecture that the difference is 
caused by the lattice corrections above T=0.6J. Indeed, ^ for the Heisenberg model is quite 
close to the value given by Eq.(|^) at T=0.6J ||^ and the deviation occurs only at higher 
temperatures. On the other hand, our data unambiguously shows that such corrections do 
not modify universal scaling functions of Eqs.(H,P). Note that were the spin waves absorbed 
by the electron-hole background in this doping range, 1/^^ and not 1/,^ would be linear in 



temperature [|13],|14|. Further, in the t-J model 1/^ and not 1/.^ is linear in temperature 



for T^J. The measurements of Ref. [|T5l also indicate linear ^~^(T, x) with nearly doping 
independent slope. 

Once it is established that the spin dynamics of doped antiferromagnets is described 
by Eq.(^), we can now address the experimental result of Imai, Slichter, and collaborators 
that at high temperatures the spin-lattice relaxation rate, 1/Ti, is nearly doping and 
temperature independent in the doping range 0-15%. From Eq.(^ one obtains: 1/Ti ~ 
/ dqXimuj-^Q'x'' {c[,uj) /uj ~ XqH"^- Indeed, we find that XqI^ and therefore 1/Ti is nearly 
doping and temperature independent in a broad doping range. Earlier, it has been shown 
that near the critical point, temperature and parameter independence of 1/Ti follows from 
the Josephson hyperscaling law; our results show that such arguments are applicable even 
for substantial doping. 

Now we turn to the Gaussian component of the spin-echo decay rate, I/T2G', which is 
given by (after Pennington and Slichter, [p!6|): 

1 0.69 2 M'q.qrp2X(q) /o^ 
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where is the static susceptibihty and Fg hyperfine formfactor. We evaluate this quantity 
by generating series directly for T/T2G and then using Pade approximants. We take J = 
1500K (in the region of experimental comparisons, calculated I/T2G is not sensitive to the 
choice of J) and the values of hyperfine couplings determined for YBa2Cu306.63 in Ref. |p!7[ . 



The motivation for doing so is that a number of experiments indicate that the hyperfine 
couplings do not change substantially from La2Cu04 to YBa2Cu306+x- Our results for I/T2G 
are plotted in Fig.^ along with the experimental data of Imai, Slichter, and collaborators 
IP for the insulator. The results are in excellent agreement with the experiments; note 
that no adjustable parameters are used. We also find good agreement at high temperatures 



(T> J/2) with the 4x4 cluster calculation for x = |jT9[. 

Recently, linear in temperature T2G above 200K has been reported by Takigawa , which 
is consistent with the QC prediction in the range of comparisons. The slope of the linear 
high temperature part of the data is larger in YBa2Cu306.63 than in La2Cu04; in our 
study, we indeed find that the slope increases with increasing doping and t/J ratio (Fig.|^, 
inset). 

To summarize, in our high temperature series expansion study we find that spin fluctu- 
ations in the t-J model exhibit quantum critical scaling behavior. Particularly, for a broad 
range of t/J ratio, doping, and temperatures (1) the numerical data for both S'(q) and x(q) 
collapse to the universal scaling function computable in the 0(N) QNLa model, (2) the char- 
acteristic energy scale for spin fluctuations, lj, is proportional to temperature and (3) 1/^ is 
nearly linear in temperature. The disagreement of the slope of 1/^ with the prediction based 
on fully renormalized T=0 value of c is likely to be caused by the lattice corrections above 
T=0.6J. One may speculate that the absence of such corrections to the scaling functions of 
Eqs.(0,^) is related to the fact that in the QNLo" model they are much less dependent on 
the cutoff and N (for a relevant discussion see Ref. [M). Our results are consistent with 



the conclusion of Refs. that in the Shraiman-Siggia model the presence of fermions 
does not alter the universality class of the z = 1 critical point separating phases with Neel 
and short range order. While the t-J model might not be quantitatively applicable for the 



doped cuprates, we believe that a broad region of the quantum critical behavior is a general 
feature of moderately doped antiferromagnets. 

In the second part of our study, we calculated the spin-echo decay rate, I/T2G, recently 
measured in a number of the high-Tc compounds. Our results for x = (i.e. in the Heisenberg 
model), obtained with no adjustable parameters are in excellent agreement with the experi- 
ments in La2Cu04 and show that the hyperfine couplings do not vary substantially from 
YBa2Cu306.63 to La2Cu04. We also find qualitative agreement with the experimental data 
on I/T2G in YBa2Cu306.63 Hi and show that for moderate doping, 1/Ti is nearly doping and 
temperature independent, in agreement with both the experimental data in La2-xSrxCu04 
P and earlier analysis based on the QNLcr model 0. 

Our results provide additional support to the conjecture of Ref. 0] that the high temper- 
ature quantum critical behavior is a common feature shared by many of the cuprate 
superconductors. More detailed study of this and related subjects will be presented in a 
subsequent publication |I0[] . 

We are grateful to T. Imai, CP. Slichter, and M. Takigawa for communicating their 
experimental results to us prior to publication, and to R.J. Birgeneau, S. Chakravarty, A.V. 
Chubukov, E. Dagotto, M.R Gelfand, L.R Gor'kov, T. Imai, D. Pines, W.O. Putikka, S. 
Sachdev and CP. Slichter for discussions. This work has been supported by the NSF Grants 
DMR89-20538 through the Materials Research Laboratory of the University of Illinois at 
Urbana- Champaign and DMR90-17361. 
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FIGURES 



FIG. 1. 

Universal scaling functions xl'^) (• ^ind solid line) and S{k) (o and dashed line). Symbols 
are obtained by collapsing our numerical data for the t-J model. Lines show analytical 
predictions for the 0(N) QNLcr model; results of Ref. indicate that there is virtually no 
dependence on N for > 3. Note that the lines are universal scaling predictions rather 
than fits. 



FIG. 2. 

^^^(T, x) for t/J=l and (o) the experimental data of Ref. [jl5[ for the undoped La2Cu04 
(we set a=l); our calculation falls within not shown experimental errorbars. 



FIG. 3. 

The spin-echo decay rate 1/T2g(T) in the Heisenberg model (i.e. for x=0) calculated 

without adjustable parameters and (•) the experimental data of Ref. for La2Cu04. Inset: 
T2g{x,T) for t/J=l. 
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